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The application of a previously developed theory of inert gas ab- 
sorption is here outlined. Interpretation of uptake curves and the 
method of obtaining tissue constants from such curves is discussed, to- 
gether with illustrations from actual experiment. A critique of previous 
analytic procedures is included. 


Introduction. In a previous paper (Smith and Morales, 1944), 
we have presented a theoretical treatment of inert gas absorption by 
the n tissues of a limb. It was there shown that the total amount of 
gas taken up to time t, ¢(¢), is given by: 


( Vi =| 
¢=C Vot—t+---+— 


Oy On 


(1) 
— (Qe + Qi eh +. +Q, omy] =C(0 -y) 


where the A’s, Q’s, and k’s are constants dependent upon certain 
physiological parameters: V, , the blood volume; FR, the delivery rate 
of blood flow; C , the delivery concentration of gas in the blood; and 
the volume V; , blood surface, S; , gas solubility, a; , and permeability, 
h;, of each tissue. It is our present concern to show how equation 
(1) and other implicates of the theory can be used in conjunction with 
experimental data. 

Fitting and Interpreting Total Uptake Curves. Let us first con- 
sider the situation in which V., R, h;, Si, a:, and V; are constant. 
The problem of fitting is then the problem of obtaining the Q; and 
the k; from the smoothed experimental data. 


* The material in this article should be construed only as the personal opinion 
of the writers and not as representing the opinion of the Navy Department of- 


ficially. 
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Regrettably enough, it cannot be said that the procedure* as we 
use it, is effortless; however, it is perfectly straightforward and al- 
ways yields results. We shall illustrate it for the important case in 
which n turns out to be 2 

1. First an asymptote is guessed, and so it is essential that the 
particular experiment be carried out long enough to permit this. 

2... From this asymptote the experimental values. of log y can 
be computed and plotted accurately on a large graph. In our experi- 
ence this plot has the general nature suggested by the heavy line of 
Figure 1. It will be found that the “straightness” of the right-hand 
portion of the curve is exceedingly sensitive (0.25%) tothe value 
chosen for the asymptote, and with a good initial guess, at most one 
trial is required to strike the proper figure. 

3. The terminal portion of the line is taken to be the graph of 
Q’e-’t, hence from this line extended and the experimental curve, it 


Log ¥ : 


F1GURE 1 


Log (¥-Qe-M) 


t 
See 2 


7 * xcept a in the fartipats but infrequent case of n= 0. 
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is possible to find y — Q’e*’*. A plot is now made of the logarithm of 
this function. It has the general nature of Figure 2. 

4. The terminal portion of this second plot is taken to be the 
graph of Q” e*"* (where Q” <0). Hence from this line extended and 
the curve, log (y — Q’e*’), it is possible to get log (py — Q’e*! — 
Qren) . When n = 2, this latter function turns out to be a straight 
line, and in fact represents Q” e*’”t (where Q” > 0). 

5. The foregoing procedure, if carried out with reasonable care, 
has yielded us results which fit the smoothed experimental curve 
throughout within the rather remarkable accuracy of better than 99%. 
Experimental data obtained by.the authors and graphed in Figure 3 


THEORY 
© EXPERIMENT 


FIGURE 3 


serves to illustrate this point. While we have not seen fit to do so, the 
precision of choice of fitting constants can then be increased by the 
standard method of alternation. Thus each component in turn is alter- 
nately computed when all the others are regarded as known. Such a 
process can be shown to converge on the best values of the constants. 
In any event the excellence of the fit, for any assumed 7, is one jus- 
tification of the theory. i 

A second method by which the theory can be tested has to do 
with the changing shape of ¢(¢) when the quantities, V,, FR, hi, Si, 
a;, and V; are changed, as few at a time as possible. How rapidly ¢ 
will reach its asymptote is governed chiefly by the &;; therefore, a 
given physiological situation will cause the steady state to be reached 
rapidly or slowly according as it increases or decreases the k;. The 
matter is schematized in Figure 4. It has been shown previously that 
increases in R, h;, Si, or a; all tend to increase the k;, while in- 
creases in V, or V; tend to diminish them. We consider cases in 
which one or the other of these conditions will prevail: 
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t 
FIGURE 4 


R is essentially cardiac delivery, and therefore should be propor- 
tional to cardiac output. Experimentally it can best be altered alone 
(to comply with the requirements of a partial derivative) by stimu- 
lating the vagus branches to the heart. 

V., as already stated, is the blood volume of the limb in ques- 
tion. The local injection of epinephrine changes V, virtually inde- 
pendently of all other quantities. 

V, changes as we compare individuals of varying “build”, say 
different degrees of obesity. These quantities and their measurement 
are discussed elsewhere. 

h, can be altered alone by the action of certain agencies, e.g. 
saponin or Ca*. 

S, can not be altered alone save possibly in development (as dis- 
tinet from growth). 

a; is altered for the case of fat as the type of fat deposited by 
the animal is changed. 

The changes enumerated above constitute possible tests for the 
theory. Granting the theory, the logic can be reversed, and a change 
in ¢(t) due to unknown causes can be analyzed. Such changes for in- 
stance are observed in fatigue or following arterial occlusion. It is 
obvious that in most of these instances more than one of the experi- 
mental parameters will be changing, consequently, in each case the 
graph must be judged specially. We are presenting such cases else- 
where. 

Absolute Constants Obtained from Experimental Curves. Atten- 
tion must be called to the important relation resulting from a pre- 
viously established (Smith and Morales, ibid) Theorem II, viz., if the 
k; are ordered, then any two successive k,;’s set the upper and lower 
bounds of the quantity, 
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for some one of the tissues. There is no particular difficulty in de- 
termining the experimental value of « for any tissue, once the nature 
of the continuous intracellular phase is known. Therefore, knowing 
either h or S/V it becomes possible to estimate the other. The quan- 
tity h of inert gases has been measured for high polymer membranes 
(Muller, 1941). The ingenious method of I. Gersh, as yet unpublished, 
yields values of S/V for various tissues. It is therefore possible to 
check and use the aforementioned relationships very readily. Such 
calculations will appear elsewhere. The ratio of the permeabilities of 
two different gases can be estimated precisely from experiments on 
the same individual using two different gases. In this application, 
suppose that a’; , h’;, S’:/V’; for the gas (’) has been found to lie be- 
tween two known k’s, and therefore have an approximate value, 
M’. If similar measurements for some other gas (”) are made, then: 


or, 


The foregoing techniques are other instances of situations where 
it is possible to glean information regarding cellular constants from 
experiments performed on the intact animal or human subject. 

Other Equations. We shall close this paper with a note on other 
equations which have been used in the study of this same phenome- 
non. 

The equation of H. W. Haggard (1924), 


i Beene ee 4B constants (2) 
o(B— 64) | 
is frankly empirical, and as such has served a useful purpose, espe- 


cially in studying the uptake of ethyl ether by men. 
The equation of A. R. Behnke, 


t 
dd —k or ¢=99(1— e-kt) (3) 
i rea " 
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was originally proposed* by Behnke as a provisional empirical equa- 
tion, and in that capacity has been very successfully applied to the 
exchange of nitrogen in human subjects and in dogs. An attempt has 
been made, however, to give theoretical interpretation to this rela- 
tionship, and such studies have been accepted widely enough to merit 
critical examination. 

It is to be noted that equation (3) is a special case of equation 
(1), for n= 0, i.e., for the case of a single absorbing phase. That the 
body cannot, with respect to any gas, be considered as a single ab- 
sorbing phase has of course been realized by all authors. However, 
it has been assumed that although the limb is a multiphase system, 
the various phases differ in their absorptive powers and therefore 
take up the gas in sequence, so that at any one stage equation (3) 
will hold, with ¢° as the asymptotic amount of gas in the tissue. 

First, there is a serious theoretical fallacy involved in this pro- 
cedure, namely that a reaction cannot depress its rival unless the 
attending A u is vastly more negative than that of the rival, and such 
is not the case in this instance. This difficulty leads to a number of 
paradoxical consequences: 

The “elimination constant” varies continuously over the range of 
time; so that there is little advantage to studying k(t) instead of 
$(t). 

The average elimination constant for any one tissue is calculated 
over a period of time when presumably only that tissue is absorbing. 
Yet when these same elimination constants are used together, the 
stages during which these tissues are absorbing are found to overlap. 

Q’s sometimes turn out to be negative numbers, and therefore 
could not possibly represent amounts of gas. 

Second, there is the important anatomical argument that the 
blood distribution of a limb is such that some blood is at all times 
the exclusive monopoly of any one type of tissue, hence it is unavoid- 
able for tissues to absorb simultaneously. When (as has been done) 
the whole body uptake is considered, the situation is even more criti- 
cal. For example, in the mesentery and in the retro-peritoneal region 
there are depots of fat which have their own vascular beds, complete- 
ly independent of the blood supply to other phases. 

Finally, it should be noted that an equation cannot reasonably 
be regarded as a theoretical interpretation until the constants which 
it involves are related to well-defined physical quantities. It is this 
limitation on the part of equation (3) that precludes the correlation 
of changes in experimental set-up with changes in “k”. 


* Personal communication. 
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The results of a simple analysis of the situation in the lungs with 
regard to the uptake of inert gases is presented. It is shown that the gas 
transfer at the lung governs the early stages of gas uptake by a body 
region so that the latter conceivably could be used as an index of lung 
function. It is shown that the so-called Fick Principle, which serves as 
the basis of most indirect cardiac output methods, is an approximation to 
a more general equation. When the latter equation is applied to a slowly 
Pere gas, it offers the possibility of determining functional lung 
surface. 


Being the spatial intermediates between the atmosphere and all 
body tissues, the lung membranes potentially govern all inert gas ex- 
changes.' It is, therefore, unavoidable that their role be included in 
such exchange theories as have been attempted (Smith and Morales, 
1944a, hereinafter referred to as SM). Furthermore, the kinetics of 
inert gas exchange across the lung membranes has been made the 
basis for certain indirect studies of circulatory function, and so, from 
this point of view, a study of these kinetics has immediate practical 
application. In appreciation of these two incentives, the following 
note attempts to formulate certain transfer problems at the lung, and 
then to incorporate the solutions into more comprehensive theories 
and applications. 

As must forever be the case in biology, the physical situation in 
the organism must be considerably simplified and schematized before 
the problem involving it can be formulated. For our present discus- 
sion the theoretical model is pictured in Figure 1. We imagine the 
entire pulmonary circulation to pass through a lung blood chamber of 
volume, V cm’. The walls of this chamber are the lung membranes 
Ce eee ee tng thar chbinisat the Navy Dopacrmant 


officially. 
1This is not to imply that only lung membranes may pass an inert gas. 
Behnke (1941) has shown that N, can be exchanged by the skin in significant 


amounts. The same author (1938) has also demonstrated that He is passed by 
the bladder wall in very slight amounts. 
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: as VASCULAR SURFACE= Si 
Ve _— 


LUNG PERMEABILITY=h, 


SYSTEMIC CIRCULATION 


FIGURE 1 


themselves. On the external side of the membranes the gas in ques- 
tion is assumed to exist at a mean constant. pressure of p dynes cm”. 
The remaining symbolism is either obvious from Figure 1 or is con- 
sistent with that previously introduced (SM). 

Gas exchanges at the lung can be divided conveniently into those 
processes which result in a net gain of gas by the body as a whole 
(saturation), and those processes and phases which result in a net 
loss (desaturation). Not only has the chief experimental concern of 
the authors been with the former process, but also the analysis of 
lung desaturation, even at the level presented here, is beset with a 
great many complications which at this point we have found it im- 
possible to simplify. For these reasons the attention of this paper is 
focused on saturation phenomena, and mention of the converse pro- 
cess will be slight. Some measure of simplicity is attained by divid- 
ing saturation itself into certain temporal stages, and this accord- 
ingly has been done. We now turn to the processes themselves. 

_.- Saturation, The time course of saturation may be pictured: as 
follows: At the instant that the subject takes the first breath of the 
gas, he builds up the alveolar pressure of the gas to the value p. We 
designate that instant as zero time. The blood at that moment in the 
lungs will then begin to flow out with a progressively increasing con- 
centration of gas. Finally, the last part of that blood will flow out. 
with a concentration of gas to be sustained for as long as p == Cone 
stant. This initial stage, equal to the length of time which blood. 
spends in the lung, we may designate as stage (a). There will now 
follow a-stage in which the arterial blood is carrying and delivering 
blood to the various absorbing regions, and during which the venous 
return does not yet carry gas. This stage (b) ends when the veins. 
HCN begin to return charged blood to the lung, and it is, therefore, 


2 erally we will Présiese the coronary return ( about 5% of the total) which 
is the first such return. 
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equal in length to the least round trip time minus the lung time. From 
this time on there is a stage (c) during which more and more re- 
gions are contributing to the venous return, ending finally when the 
region with the longest round trip circulation time just begins to con- 
tribute to the venous return. From this instant on, for as long as the 
experiment is in progress, all regions are contributing to the ve- 
nous return. We call this last stage (d). During each stage we will 
want to know: (1) the concentration x of the gas in the blood leav- 
ing the lung, (2) the concentration y of the gas in the blood arriving 
at the lung, (3) the uptake © and the rate of uptake ©’ of gas by the 
body as a whole during this stage. The time relations of the latter 
variables (3) are schematically depicted in Figure 2. The nature of 
the curves will be discussed in the analysis. 


RATE OF UPTAKE 3° 


UPTAKE, & 


(a) (b) (c) (4) 
FIGURE 2 


Stage (a). Letting x denote the concentration of the gas in the 
pulmonary veins at time ¢, we see that x will grow in time according 
to the same law as it would have grown had the blood in V been stag- 
nant, the only difference between the two cases being that in the for- 
mer instance x will never attain its asymptotic value, but will be cut 
short and kept at a reduced constant level. This follows easily from 
the equation for solution of the gas into the moving blood. Due to the 
fact that in the lung the blood flows rapidly through exceedingly fine 
vessels, we may assume that diffusion within the blood in a plane 
transverse to the axis of the vessel takes place instantaneously, while 
no diffusion takes place along the axis of the vessel at all.. To this 
approximation, with k as the Henry’s Law solubility constant,’ 


: 3 Thi tion is actually derived by imagining the blood. stream. split.into 
finden ements of elas thatau through V at a rate such that FR cms® will 
S aditough V per second.. The concentration of gas attained in each of these is 
hee calculated, and from this, the concentration of gas in. the emerging, blood. 
Inherent in the latter step is the assumption that the ratio of surface to volume 
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V dx/dt=hS (kp — x). (1) 


Under the initial conditions, 7 = 0 when t = 0, we integrate equation 
(1) to give 
2=kp(1 — e*8/"), forOStSV/R. (2) 


Thus x starts from zero and rises to 
2=kp(1 —e"*/2), for V/R St, (3) 


at which value it will stay so long as p and R remain constant. 
During stage (a) the rate of uptake of gas is 


®' (g) = Rx = Rkp(1 — e’*’"), (4) 
and the uptake itself is 


t t 
©) = f e'wat= { Rkp (1 — e*8*/") dt 
0 0 


(5) 
= Rip |t V(t emery /ns| 


Stage (b). Passing into the next stage we see that x remains 
at the value indicated by equation (3), whence 


3, = Rx= Rkp(1 — e"*/*), (6) 
and 


by) i Rkp(1— ewe) |( t -— ) (7) 
Much importance attaches to stage (b) for two reasons. First, 
it represents the stage at which tissues commence to absorb gas, and 
second, it is the stage during which cardiac outputs are usually meas- 
ured. We will digress to follow these two considerations. 
It will be seen that the delivery concentration, C‘”, at a distant 
region, ”, will vary in the following way: 


Co=0, fordOstst,.™, 
Ce kp (1 — e*/7), fort,9 StSt,0+V/R, (8) 
Co = kp(1 — e*/®) = xo), fort.” +V/RStES ow, 


where t,‘” is the “arterial time’, or time taken for the blood to flow 
from the lung to the 7-th absorbing region. Since in the normal mam- 


_—__oooO::::- >: ere—a SS  ——— ee 


of an element, AS/AYV, is the same as S/V. To the extent that all blood vessels 
are roughly cylindrical save at junctions, this assumption would seem justified. 
We are also tacitly assuming that the flow of blood is frictionless since we assert 
that a local concentration can be sustained despite the streaming of the blood. 
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mal both ¢, and V/R are a matter of seconds (see Kvale and Allen, 
1939, and Schmidt, 1941), as contrasted with a gas saturation time of 
over three hours (Smith and Morales, 1944b), it is not unreasonable to 
suppose that C™ builds up instantaneously to x(), although this is 
of course not strictly the case. However, in certain situations, of 
which a most important one is pulmonary edema, or again in shock, 
both ¢,°9 and V/R may be significantly lengthened due to an im- 
paired circulation. In such cases C“” can no longer be considered con- 
stant in the early stages of absorption. These facts have a certain 
bearing on the gas exchange of a limb, and more specifically, on the 
shape of the limb absorption curve. Thus, the differential equations 
governing absorption in the 7-th region are (SM) 


Vo da” /d§é= RO CM — RO a, — > hiSi (to — 25), 


=1 


: c= BOLE) © Rie ae ee omer e ter eo) 
Vy” dz; /dt= hS,™ (29°? — a," ), 
and the total uptake up to time t is 
fp = Vo 4.69 + V7, +.--+ Vi a, (10) 


while 
dp =R” (cm = GoiO) 


The early variation in C’ can be taken into account exactly by solv- 
ing equation (9) in stages as prescribed by equation (8). While that 
is not a formidable problem, its explicit solution is fairly tedious, and 
for the time being we have not the experimental equipment rapid 
enough to verify the result on normal animals. However, the general 
quantitative description of this phenomenon is easily perceived with 
the aid of equation (10). For when t = 0, x‘ =0, so that 4” (0) 
is equal to R™ C™. In the approximation, C’ = constant > 0, so 
that (0) > 0; in the exact solution, C‘ (0) = 0, so that ¢’ (0) 
= 0, Since in either case 4’ (co) =0, while 4“ (4) > 0, (0<t< 0), 
it follows from Rolle’s theorem that the exact solution for ¢'” (t) 
passes through at least one maximum, and, therefore, that go (t) 
passes through at least one inflection. These relationships are depicted 
schematically in Figure 8. In a number of uptake curves we* have 
found it to be the case that a hint of an inflection occurred within the 
first minute of absorption. It is to be inferred from the above that 

i i i i iscuss this 
eo ae Gi yceaity op Califortins, Fla 'tod bad observed 
this phenomenon in several carefully determined curves and had successfully 
demonstrated their reproducibility some years ago, but his work has unfortunate- 


i i is 1 i ly test of 
emained unpublished. The presence of this inflection as a very ear. 
He ayahinetion remains a very interesting possibility. 
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such an inflection will be accentuated and will appear later whenever 
either or both pulmonary and systematic circulations are abnormally 
slow. 


C or > 


oO | minute time — 90 minutes 


FIGURE 3 


As remarked above, the most useful methods of determining car- 
diac output, namely those methods using an inert “foreign” gas such 
as acetylene, are put in operation in stage (b). In essence, these 
methods all depend upon the determination of #’ from measurement 
on the rate at which the inert gas disappears from the spirometer, 
and simultaneous calculation of ®’ as = Rx. 

A convenient method of doing this for continuously determinate gases has 
been communicated to us by Dr. H. B. Jones of the University of California: 

Since, from a finite breathing mixture of oxygen and inert gas the respira- 
tion will progressively remove more oxygen than inert gas, the concentration of 
inert gas relative to the concentration of oxygen would normally increase in time. 
However, it is possible to determine experimentally how much oxygen per unit 


time must be added to the respirometer in order that the ratio of the concentra- 
tions of the two gases will remain constant. : 


If K° = initial volume of inert gas in the spirometer, 
va = initial volume of oxygen in the spirometer, 


ae = oxygen consumption per unit time, 
B = oxygen addition per unit time, 


then it can be easily shown that 


Ko +t 
pe ya he 

yo? 

Os 0 


It re been found experimentally that during tha atape (b), B is a constant; 
consequently, for this stage, we may write 


i Gps “ake 
*=(5 a )exor, 
v 
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where the coefficient of ¢ is fully determined experimentally. Combining with 


equation (18): 
Rea /o We 
={( Soe ) p. 
Oa 


In the same subject the quantities p and k can presumably be kept constant, and 
in that event, it follows that (K°q/v? ) — £8 is a measure of cardiac output. 


Jones has used this measure clinically with encouraging results. 


Then, of course, 
R=0./7. (11) 


Substituting from equation (6), we see that as long as we are truly 
im stage (b), 
R= 4'/kp(1 — e"8/*) = cardiac output . (12) 


Knowing the various constants and determining ©’ experimentally, 
Wwe may solve equation (12) which is seen to be transcendental in R. 
However, if we furthermore assume that e-"’/® ig negligible, then 
equation (12) reduces to 

R=®6'/kp. (18) 


Actual methods make both assumptions (12) and (13), and so it be- 
comes necessary to study these proposals and their verification. 

It is be noted at the outset that at least for C,H, Grollman (1932) 
has very conclusively proved‘ the truth of assumption (13) by com- 
paring F values so calculated with R values obtained by the so-called 
direct heart puncture method. And it is a further suggestive fact 
that the errors introduced by the failure of the two assumptions tend 
to cancel each other (both numerator and denominator of (13) are 
smaller than they should be). 

We may answer the question of neglecting the exponential by 
recourse to data in the literature. The basal level of R (Grollman, 
ibid.) is about 66 cms* sec. The quantity S has been variously esti- 
mated to be 5 to 9 X 10° cms”, of which measurements perhaps Will- 
son’s (1922) is the best, namely, 7 < 10° cms”. Thus, even if h were 
of the order of 1 X 10~ cm sec", the error introduced by omitting the 
exponential would be under 5%. Actually, it is to be expected that 
h will be greater, because histological examination of the alveolar 
wall has failed to show the presence of a continuous plasma membrane 
Fee ee eee sioettaey (LOAD) Waa qemered 
considerable evidence to show that the lung may contain variable blood depots 
whose contents are not at one time a part of the pulmonary flow. It is evident 
that if some of the gas inspired (as measured by the respirometer loss) over the 
trial period were absorbed by the depots, then the A-V oxygen difference divided 


into the respirometer loss would give values of R which are too high. This error 
would be inherent in all gas methods. 
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between the gas of the lumen and the blood (Maximow and Bloom, 
1940). Translating the Kroghs’ “diffusion constants” into absolute 
permeabilities, we have obtained® 


Boyer 7.00 X 10-° cm sec?; heo, = 8.0 X 10+ cm Ssec?, 


while for oxygen and carbon dioxide respectively, the e”’/* values 
are .0000 and .0408. Thus, these two gases would be satisfactory in 
cardiac output measurement if they could be regarded as inert gases. 
In any event, it is clear that a determination of the h of the gas is 
prerequisite to deciding in favor of either equation (12) or (13). The 
latter equation, (13), has also another experimental implication. If 
one chooses a gas for which h is small, so that the exponential is not 
negligible, we obtain instead of expression (13) the following: 


Rkp — & 
Rkp 


The quantity ®’ for the gas is readily obtainable, while R can be de- 
termined either directly or by the acetylene method. From the above 
formula we may, therefore, compute hS , or if h can be independently 
obtained (as it can, although in a more difficult measurement), S, the 
functional (as distinct from morphological) surface area of the lung 
can be determined. This is a quantity not otherwise easily obtained, 
and has considerable importance in metabolic theory. 

An interesting special case of gas absorption is suggested by M. Krogh 


(1914). For gases which are rapidly and irreversibly combined with some fixed 
constituent of the blood (e.g., hemoglobin), we have in stage (a) 


hS = R log 


emo + orig eS : (i) 
V 
where x° is the original concentration of gas in the blood. In stage (b) we reach 
hSkp 
(00) =a + ——. (ii) 
The rate of uptake during stage (b) is then 
&'’ —hSkp, (iii) 
while the amount lost from the respirometer is 
& = (hSkp)t. (iv) 
On the other hand, the total amount present in the body is 
(Original amount) + ®= va + (ASkp)t, (v) 


5 See M. Krogh, 1914 and A. Krogh, 1919. It is interesting to note that the 
measure of permeability used by these authors is stated in terms of a pressure 
gradient across the membrane, while absolute permeabilities are written in terms 
of concentration (in solution) gradients. Due to its high water solubility, the h 
value of CQ, is lower than that for oxygen, contrary to the widely held view. In 
our opinion, the latter is a more natural measure, since the gases certainly per- 
‘meate only in solution. 
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where » is the total body blood volume. Relation (ii) can be made to be an equa- 


tion for cardiac output, and the slope and intercept respectively of (v) can yield 
pand hS. 


For pointing out a possible fallacy in this equation we are indebted to Lieu- 
tenant N. Pace, H-V(S), USNR. If the equilibrium between CO and HbCO so 
overwhelmingly favors the latter that all the CO available is consumed without 
exhausting the carrying power of the blood, then R could not be inferred from 
measurements on the disappearance of gas from the external circuit. This con- 
dition, however, seems unlikely in view of the fact that there is considerably more 
CO in the expired air than can be accounted for by “dead space”. 


Stage (c). With the onset of this stage, y , the concentration of 
gas in the pulmonary artery is no longer zero. The concentration in 
the pulmonary vein, x(.), of course remains constant. Thus, we have 


&',— R[x(o) — y], (14) 
and 
t 
So—- he eo b — t.)) — 6h ydt . (15) 
t-() 
Therefore, the chief problem lies in the determination of y. Let us 
proceed by arranging the N absorbing regions, r—1,2,.---N, in 
order of increasing (round-trip) circulation times. ¢,‘”, i.e., 


PO FOP) pi PO), 


Referring to zero time, it is clear that stage (c) commences at 
t=t,.™, and ends att =t,™. Blood arriving at any instant, t, from 
the 7-th region brings a contribution from that region as of 1t,‘” sec- 
onds back, at which time that region has been absorbing for ¢ — ¢,‘” 
seconds. Here we use t,‘” to represent the “venous time”, or blood 
flow time from the r-th region to the lung. If the blood in the vein 
draining the 7-th region flows out at a rate of FR‘ cm* sec, and is 
charged to a concentration of x)‘ gm cm“, then its contribution to 
the pulmonary artery concentration (when such blood arrives there) 
will be (R“/R) x,‘ gm cm-, and the amount which has arrived there 
from the r-th region will up to time ¢ be 


if Boa? dt. 
t 


XG) 


Substituting from equation (15), SM and integrating, this becomes 


Ra (0) Cl aN. Ao;” Eki (t-te) ) Qt 
t.(n) j=0 
(16) 
jen A, (7) . 
= R”2x(0) | (¢—7,) — > [1 = epeirete I} ; 


PCH) 
j=0 ke, 
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And combining expression (16) with (15), we see that for Eg aad 
<a tr, 
BD (¢) = Rx(o) (t a Eo()) 
(17) 


j (r) 
jen A. tr CL — erie t , 


St pin = bee He 

E RiP a() \(E= #4) — BG 
Thus, it is seen that in stage (c). the rate of uptake begins to decrease 
in a fashion which can be described by sums of exponentials. If the 
body could be regarded as consisting of only a few functional regions 
with radically different t,’s, then strictly speaking ®. would be a dis- 
continuous curve during this stage; however, since there are many re- . 
gions with not too different circulation times, the discontinuities are 
likely to be completely masked.* It may be noted that the greater the 
kj“ are, the more rapidly will ® vary during stage (c), and the great- 
er will be the error introduced into a cardiac output measurement 
overlapping this time stage. 

Stage (d) is easily disposed of by putting the upper limit of the 
summation in equation (17) as7—=N. Theoretically as well as prac- 
tically, ®,,) is a continuous curve. Combining the equations of stage 
(d) with previous work (SM) we see that the stage ends with 


Xo” (00 ) = 200), 
Li' (co) =X(00)/ai,t#0, 
G'(o) =0, (18) 


®(00) = 2 (00) (9+ S SVi/a). 


f=1 I=1 


The state described by equation (18) is theoretically never 
reached, for ® is asymptotic to 6(0«). However, within a few hours 
of gas uptake ® is experimentally indistinguishable from ®(o), and 
for practical purposes we speak of the body as being “saturated” with 
gas. 


This completes our analysis of the saturation process. 


Desaturation. If while equation (18) holds, p is suddenly re- 
duced from » = constant to p = 0, the entirety of the resulting pro- 
cesses is termed desaturation. In principle, this process is governed 
by the same general equations as its converse. 

By solving equation (1) with p = 0, we have immediately the 
relationships between the concentration of gas in the pulmonary ar- 


_ © Another agency which in actual practice ensures this occurrence is the mix- 
ing that results from the velocity profile within blood vessels. This has been re- 
ferred to previously. 
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tery and the pulmonary vein: 
% = ye-s/V)t 0<t<V/R, 
L=YyewsR . V/R<t 


The tissue discharge is still governed by equation (9) with ap- 
propriate boundary conditions. In particular, if C =0, then 


x” = x (co) (Ap et ar Ane™ + see + Aone **), (20) 


As usual, the larger the k,’s the faster will be the variations in z. 

Finally, dilution at the arterial end of the lung takes place ac- 
cording to the same law. Denoting by x,‘ (¢ — t,‘) the value of 
Xo‘ ty” seconds ago, then at time ¢ the concentration of gas in the 
pulmonary artery is 


(19) 


1 = 
R ay R24” (¢ = t,("). (21) 
r=1 


What introduces a hopeless cumbersomeness into the analysis of 
desaturation is the fact that different regions have different circula- 
tion times; hence, all regions are always at different stages of dis- 
charge, and remain so for all time. In equation (19) the exponential 
appears alone, so that unlike the situation in equation (12), it cannot 
be disregarded (there is ample experimental evidence to support this 
assertion). More so than in the case of saturation, the existence of 
portal circulations precludes any such assumption as t,“) > €, im- 
plying t,@ > ¢,. Thus, as a whole, any simplifications would be dis- 
tinctly out of accord with experimental facts, and we have not seen fit 
to press the analysis further. 

Notwithstanding the difficulties of interpretation, the desatura- 
tion phenomenon is clinically important, and has received consider- 
able experimental attention. Interesting contributions have been made 
by Behnke (1935) on N., and by Underwood and Diaz (1941), who 
injected radon into the saphenous vein and by its radioactivity meas- 
ured its appearance in the expired air. The latter authors present this 
data in the suggestive fashion of the plot, log (—®) vs. ¢. The ex- 
periments were carried out over four minutes, yielding an average 
point every minute. It is quite clear that in every graph shown, the 
last three points lie on a good straight line, while the first is always 
displaced upward. This would seem easily explainable in a qualitative 
way: On the first lung passage and before the next (roughly one min- 
ute), the loss is governed by expression (19) with only slight dilution. 


7 Note that equation (19), like equation (13), is an equation for functional 
surface area of the lungs. 
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On subsequent passages the main loss is still governed by (19), but 
with lowered concentration (due to complete mixing). Accessory loss- 
es to tissues are of course also present, and governed by expression 
(9). The value given for the first slope is about 0.69, not too different 
from, and lower than the hS/R computed above, as is to be expected. 
The subsequent stages involving tissue absorption are of course of a 
lower and uninterpretable slope. 
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ON THE FORM AND STRENGTH OF TREES: PART II 
THE PRIMARY BRANCHES 


I. OPATOWSKI 


ARMOUR RESEARCH FOUNDATION, ILLINOIS INSTITUTE OF 
TECHNOLOGY, CHICAGO 


_ The main purpose of this note is to relate the shape of the branch 
ie Se eee and to explain by theory why a tree has a large number 
of branches. 


1. The branch as a cantilever. If the weight of the secondary 
branches is taken into account by a suitable increase of the specific 
weight of the primary branch, a reasonable model of the primary 
branch is a horizontal cantilever of conical shape. (A model of a 
branch of maximum strength will be discussed in Part III.) If 7’ is 
the radius of the base of the cone, L the length of its axis, p’ the ratio 
of the weight of the complete secondary branches to the primary 
branch, m™’ the ratio of the weight of the primary branch in green 
state to its weight when oven dry, the maximum bending stress, 
which occurs at the base of the cone, is (Timoshenko, 1940, pp. 88- 
93): 

m'(1 + p') GL?/(30007’) (1) 


where G is the specific gravity of oven dry wood. The units are 
throughout the paper kg and cm unless otherwise stated. Putting 
fraction (1) equal to o,/f', where f’ is a safety factor and o, the 
maximum possible bending stress which is given by equation (3) of 


Part I, we get: 
m' (1 + p') f'L? = 3711000 G4 7’. (2) 


Assuming that the safety factor is the same for all branches of the 
tree, we see that for a given species the length of the branch is pro- 
portional to the square root of the radius at the base. For the species 
grown in this country G varies between about 0.3 and 0.7 (U.S. De- 
partment of Agriculture, 1940, pp. 50-53), so that approximately 
G)/4 = 0.75 to 0.91. Taking for instance m’ = 2.5, p = 0.6, f' = 84, 
G = 0.5 so that G14 = 0.84, we get for 7’ = I cm, L = 1 meter. 

In applying the theory to an actual branch it may be convenient, 
because of the irregular shape of the branch, to eliminate the length L 
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in relation (2) by means of the total weight of the branch. Therefore, 
calling B the weight of the primary branch in the state as it is ob- 
tained after being cut from the tree, that is, with all its secondary 
branches and leaves, we have 


B=am' (1+ ')Gr? L/3000, (3) 
and from expression (2) we get: 
f B?=4.07 m'(1 + p)G*r*. (4) 


Or, if P is the weight of the primary branch alone in the green state, 
that is, without secondary branches and leaves, we have: 


f(ltp)P?=407T mG’ r®. } (5) 


2. The shape and the strength. We eliminate the radius 7 in 
(4) by means of the lateral surface S of the primary branch. As 
Rashevsky (1943, p. 44) pointed out, S has a biological meaning, be- 
cause in a first approximation one may assume S proportional to the 
area of the leaves and therefore proportional to the amount of sub- 
stances exchanged by the branch with the surrounding air. Since 
for a cone of the size of a branch 1 + (7’/L)? = 1 we have, 


Sar Ls (6) 


The elimination of 7’ in equation (4) by means of relations (38) and 
(6) gives: 


10° CB’ = f'[m'(1 + p’) ]* G'’* S® , where C = 989 . (7) 


We attempt now to explain some shape characteristics of the tree 
by considerations of mechanical strength. To do this, we will have 
to compare with each other some shapes of branches and branch ar- 
rangements which are geometrically different but biologically identi- 
cal. A reasonable mathematical condition for a biological identity is 
the constancy of G, m’, p’, B, and S, which will be assumed through- 
out the following discussion. It must be noted, however, that actually 
the weight B consists of two different parts: B, = the weight of 
leaves and sapwood which have mainly a biological function, and B,, 
weight of the heartwood whose main function is to provide the tree 
with the necessary mechanical strength. The ratio B;/B, might be 
called the biological efficiency of the branch or of the tree. A mathe- 
matical condition for a biological identity would be B, = const. rather 
than B = const. But if we include into the concept of biological iden- 
tity the condition of equal biological apa By = = const. ete Ki 
substituted by B = const. 
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., We discuss now the following questions: 

_ (1) Why does the cross-section of a branch decrease towards 
the free end? For a conical shape of a branch the safety factor jf’ is 
given by equation (7). For a branch in the shape of a circular cylin- 
der we find, if we neglect the area at the free end in calculating S, 
the same equation (7) with C =.195. We see, therefore, that the 
conical shape makes the branch five times stronger. 

(II) Why does a tree have many branches? Consider a tree 
with » conical branches and assume that the safety factor f’ is the 
same for all branches. Call 7’; the radius of the base and L; the 
length of the i-th branch. Instead of equations (2), (3), and (6), we 
have now the following equations: 


r,=m' (1 + p') f'L,?/ (8711000 G) , (8) 
B=am'(1+ p)GS% 1',? L;/3000 , (9) 
S =n be Ts Li r) (10) 


where B is now the total weight and S the total lateral surface of all 
branches. If we eliminate 7’; in equations (9) and (10) by means of 
equation (8) and calculate B?/S®, we get equation (7) with C — 989K , 
where 


K= (34. L8)9/(3%, L9)?. 


K is > 1 for n > 1, whatever are the values of the L;’s. This is an 
immediate consequence of Jensen’s inequality (Hardy-Littlewood- 
Polya, 1934, p. 28). Therefore, independently of the length of branch- 
es, a tree with n branches is K times stronger than a biologically equi- 
valent one-branch-tree. To get an idea about the order of magnitude 
of K , we put L; =7L,. Then (Pascal, 1900, p. 235): 


K = (27/16) (n? + n)+/ (2 + 2n—1)*. 
Several values of K are given in the following table: 


nw 1-2 3 4 5 6 7 8 9 10 25 
K=1 16 2.8 46 6.7 92 121 15.5 19.3 23.5 137 
forn—> ©,K-— (27/128) n?. 


In conclusion, one may say that the tree shows some kind of ten- 
- dency to shape itself in such a manner as to increase its own strength. 
However, this may be a question of viewpoint only. In fact, if one 
compares a one-branch-tree with a n-branch-tree on the assumption 
that the strength is in both cases the same (f’ = const.), one finds 
[ef. equation (7)] that S is proportional to K’”*, so that the large 
number of branches that a tree has may be interpreted also as a ten- 
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dency of the tree towards a higher biological efficiency. The increase 
of S is, however, much smaller than the increase of f’ in the previous 
discussion, for instance, for n = 25 we have K’’* = 2.7. 
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THE TWO-FACTOR THEORY OF NERVOUS EXCITATION 
WITH NON-NORMAL ACCOMMODATION 
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: By introducing a parameter which vanishes when the accommodation 
is normal the modifications in theoretical predictions when accommoda- 
tion is non-normal are discussed, particularly in the case of stimulation 
by alternating currents. It is found that non-normality always has the 
effect of lowering the optimal amplitude of the alternating current, 
measured in units of the observed rheobase. 


Quantitative comparison of experimental results with predictions 
made by the Rashevsy-Hill two-factor theory of nervous excitation is 
always made with the assumption of normal accommodation, and no 
one seems to have considered in detail what is to be expected in case 
the accommodation is non-normal. It is true that the equations be- 
come quite complicated then, and the agreement is very good in a 
great many instances even with this restriction. On the other hand, 
there are deviations, especially with stimulation by alternating cur- 
rents (Hill, 1936; Hill, Katz, and Solandt, 1936; Coppée, 1986) where 
the optimal amplitude of the stimulating current is rather less than 
predicted (Parrack, 1940), and it seems worth-while to inquire 
whether this deviation—and possibly others—can be attributed to non- 
normality of accommodation. 

We shall suppose that the deviation from normality is never 
great, and we shall introduce a parameter for measuring this devia- 
tion and calculate the corrections of the first degree in this parameter 
in the case of stimulation by alternating currents. It is convenient to 
use Rashevsky’s form of the equations and to introduce the Haupt- 
nutzzeit as the unit of time. If we write y, and y, instead of « andj, 
this, together with a suitable choice of units for measuring the y;, 
enables us to write these equations in a normalized form (Household- 


er, 1944) ; 


yi tu(yi— yo) HeS, (1) 
with primes denoting the derivative and S representing the stimulus. 
If we write 

log (ti/t2) =A= fa — be +7, (2) 
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then 

f= (A— »)/(1— e&), 

fa = (A— »)/(e— 1). 


With » = 0 we should have normal accommodation. In general 4 > 0. 
The cathodal rheobase R, (Hill’s “observed rheobase’’, I,), is 
found to be given by 


HO Bo 1 1 eta ets 
Marnier pee (4) 
R, fe fh Me la 


(3) 


while that which Hill calls the “true rheobase”, J, , the rheobase in 
the absence of any accommodation, is given by 


(Yo? — Yr) /Lo = CM/ th « (5) 
In the case of stimulation by alternating currents 
S=TIsinat, (6) 


and in the steady state 
yi — YX =A. Tsin(wt — yi), 
A; = e4 (w? + i?) 7, (7) 
sin y; = @(w? + y;7)-1/?, 
COS yi = yi (w? + yi?) 7”. 
The difference 
a= vey ge Ye) 1 /L 
is expressible as a periodic function whose amplitude A is given by 
__ w? (et — em)? + (yy ef — pp eM)? 
(w? + uy”) (@? + tu?) 
and for minimal stimulation, 
AI=y,.° — y,°. (9) 


Since A is a function of w, the least amplitude J necessary for 
stimulation is that value [,, obtained from equation (9) with A maxi- 
mal. If we set 


A? , (8) 


A= (e — e)/B, (10) 
the maximal A is given bythe minimal B . 
We now set 
. F = (4, e¥ — pp em) / (eH — et), 
(11) 
M;?= Mi? = F?, Q2= w+ F2, 
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which enables us to write 


B? = (Q? + M,?) (Q? + M,*) Q? (12) 
and obtain easily the minimum 
B=Bn,=M,+M.,, Q,?2=M,M,. (13) 


We note that F vanishes with », that M;? differs from ,;2, and 
hence, M; from yu; only by terms of the second and higher orders in 
y, and that Q differs from w to the same extent. With » very small, 
this difference could become significant, but not otherwise, and we 
may regard the optimal frequency as independent of the normality of 
accommodation. Likewise, the symmetry of the graph of B (or of I) 
against the logarithm of w (Hill, 1936) is unaffected except for small 
values of » , and also the rectilinearity of the plot of B? (or J?) against 
w* for large w. On the other hand, the plot of B? (or J?) against 
w? on the low-frequency side of the optimum should deviate from 
rectilinearity for sufficiently small values of w. 

We shall see, however, that non-normality of accommodation tends 
to diminish J,, where y > 0. We note first that Hill computes the ratio 
of J to the “true rheobase’, J,. Since, with sufficient approximation, 
M; =, we find from equations (5), (9), (10), and (2) that 


De ee rte jen) / CL ie te 2?) = (LT 6) CLG es UE) (14) 


from which it is clear that this ratio always exceeds unity and in- 
creases with w. However, if we use equation (4) instead of (5), we 
find, after making an expansion, 


See 


Xr 
In/Re= (1 + e+) e-ri(e*=4) 2 a [i-=( Te) ewe ay] y.(15) 
Now for positive 4, the term independent of » increases monotonically 
from 2/e when 4 = 0 to 1 when i = ©. Since this is the case, it fol- 
lows further that the coefficient of » is always negative for all values 
of 4 whence, for a given ratio 1/1, the less nearly normal the ac- 
commodation, the lower will be the amplitude of the stimulating al- 
ternating current, at optimal frequency, in units of the observed rheo- 
base. Coppée (1936), for example, reports ratios of from 70% to 
80%, and while Hill, Katz, and Solandt (1936) regard the experi- 
mental procedure as faulty, such values may, as we see, be regarded 
as an indication of non-normality of accommodation, especially in 
view of the fact that these investigators themselves find ratios which 


are below their expectations. 
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In Figure 1 three sets of Coppée’s (1936) experimental data are 
compared with theoretical predictions, and the values of 4 and » are 
indicated. In the case of Coppée’s fourth curve, no satisfactory fit 
was obtainable, even with corrections of higher order in », so that 
evidently other factors are involved. Nevertheless, it is clear that 
low experimental values of the optimal amplitude can be quite well 
accounted for by the two-factor theory, although, of course, devia- 
tions in other respects may occur as in the case of Coppée’s fourth 
‘curve. 

In making this comparison the experimental values were meas- 
ured on the graph as published, since tables were not given, and the 
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plot was made doubly logarithmic for convenience of fitting. In this 
way it is possible to draw a set of theoretical curves with different 
values of 2, and by horizontal and vertical translations choose that 
one which fits most nearly. This determines 4. The height of the 
minimum then determines », and its abscissa could be used to deter- 
mine the time-unit. These computations were made by Mr. Manuel 
Schwartz. 


This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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An expression for the intensity of central excitation cortesponding 
to the perception of an angle is derived. A general expression for the 
aesthetic value of certain types of polygonal patterns is derived there- 
from and compared with experimental values obtained by the rank order 
method. The theory is found to agree with the experiment within a limit 
of about 10%. 


In a previous publication (Rashevsky, 1942, p. 179) we have dis- 
cussed the relation of the intensity S of the peripheral stimulus cor- 
responding to the perception of an angle ¢, to the angles 6, and 6 
which characterize the position of the sides of ¢ with respect to the 
vertical. Consider the simpler case that the stimulation is the same 
every time that the vertical (or horizontal) component increases or 
decreases by the same amount, and that the coefficient of proportion- 
ality a’ is the same for the vertical and Cube component. Then 
we find: se 


Resin in 6h scons 6 — cos 4, (ig, we! UL) 


In another publication (Rashevsky and Brown, 1944), we have 
shown that the central excitation increases as the square of the peri- 
pheral intensity S. Hence, with a as another coefficient, we find for 
the intensity of excitation of the center corresponding to the percep- 
tion of an angle ¢ the expression: 


; C49 = a(|sin 62 — sin 6,| + |cos Ba SOR Ea ate (2) 


it has been suggested eaceuky: 1942, ‘pp. “189- 190) that rela- 
tions such as equation (2) may be checked experimentally, by measur- 
ing. aesthetic values of polygons in which lengths. and. positions of 
straight lines in the plane do not vary, but angles do.. This is obtained 
by permuting the sides of a polygon, leaving, however, the direction 
and length of each side, invariant. Thus, starting with, the 16-sided 
polygon shown in Figure la, we may obtain by permutation. of the 
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a b c 


FIGURE 1 


sides such polygons as shown on Figures 1b and Ic, and a large num- 
ber of other widely different shapes. For all such polygons the num- 
ber and the intensities of excitation of all L-centers, V-centers, and 
H-centers (Rashevsky, 1942) are the same. However, the symmetry 
relations may change from case to case. Some of the polygons have a 
central symmetry, like Figure la. Some have a vertical and horizontal 
axis of symmetry, but no central symmetry. Some have only a verti- 
cal axis of symmetry. Some may have no symmetries whatsoever. A 
total of 77 polygons obtained by permuting the sides of Figure la 
were used in an experimental study of their aesthetic values by the 
rank order method. All 77 polygons had a vertical axis of symmetry, 
but not all had a horizontal. 

Consider the general expression, derived previously [Rashevsky, 
1942, equation (28) ], for the total intensity # of central excitation, 
corresponding to a pattern of n elements. With the same notations: as 
before, we have: 


E seed Nie, —* (1) Ser 0y, (3) 
with 
n= Ni. (4) 


The n elements constituting a pattern are divided into different 
groups such as the A-centers, L-centers, V-centers, etc. It is plausible 
to assume that while the inhibitory coefficients b; will be different for 
different groups, they will be the same for the centers of the same 
group. We shall therefore denote by 7m; , bx; , ex; the quantity refer- 
ring to the i-th neuron center of the k-th group. (Thus k = 1 may 
stand for the group of A-centers, k = 2, for the group of L-centers, 
etc.) Then, assuming that b;; = 6; is constant for a constant k, we 
have from equation (3): 
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=> [1 ma (a ae 1) b;.] > Niki » (5) 
k 4 ; ‘ 5 : : 
n=DTDd Mi. (6) 
Kk 4 


If we consider a set of polygons in which all symmetries are the 
same, and only the angles vary, then n is constant throughout the set. 
The sums = Mi €i are also constant; all except the one referring to 


the angles. Denoting now by b, the value of b; for the groups of 
A-centers, and correspondingly putting for that SYOUP Ni Cxi = Na Cai » 
we now find, with C as a constant: 


BH=C+ [1— (n—-1) ba] & Ms Cai- (7) 


But m4 is known, being the total number of angles, and é,; is given 
by equation (2), where ¢ is the value of the i-th angle. For every 
polygon the quantity > n, e4; can thus be calculated, and the observed 


aesthetic value should be a linear function of that quantity within a 
set of polygons with same symmetry relations. For sets with different 
symmetry relations, the values of n will be different, because of dif- 
ferent numbers of symmetry centers. Therefore, for different sets 
the slope of the straight line should be different, and the value of C, 
which is also a function of 7, should vary also from set to set. From 
the slopes and intercepts of the straight lines for sets with different 
values of n it is possible to determine the values of the different b;’s. 
The experimental procedure was the same as that used in a pre- 
vious study (Rashevsky and Brown, 1944), except that all the poly- 
gons were printed on square cards, 12 X 12 cm. The length of each 
side was 1.2 cm. A total of 109 subjects was used. The linear relation 
(7) was thus confirmed experimentally. Within the accuracy of the 
experiment it was found that the slopes and the intercepts of the 
straight lines were the same for all groups, except the ones with cen- 
tral symmetry. In the latter the slope is the same, but the value of C 
is greater. This means a very small value of b,, and a large value of 
e; for the center of perception of central symmetry. In order to 
obtain the values of b, from such experiments, it will therefore be 
necessary to use polygons with many more sides so as to increase 7” . 
Due to the above simplification, caused by the smallness of Bay 


we may use the expression 
E=C, + Cz + Cs 3 Ma Cai, (8) 


where C, = 0 except for polygons with central symmetry. The quan- 


VISUAL AESTHETICS 


166 


ux 


MW, 


@ Gand 


(h24+908 
0et oz Ol oot 06 08 04 


% Wand 
{1A 1A An Al 


"I 1 


CASHOLSELIAS HET ISO LIAS HPETIAGIEMUOSEBLIASPEZTLAELIASPETILBAYOSHVEZTLBLIOSHPETIAOPETIB HZIBLIS HETIL 


LN3IWI8IdX2 @——8 -——® 
ABOIH1 O-- -O - © 


N. RASHEVSKY AND VIRGINIA BROWN 167 


tity S ma ei is calculated from equation (2), and contains still the 
factor a. 

In Figure 2, the full line connects the experimentally observed 
aesthetic values; the broken line, those calculated from equation (8), 
with C, = 1.9; C, = 2; aC; = 0.11. The average difference between 
the observed and calculated value is 10.5%. The most frequently oc- 
curring difference is below 5% 

A remark must be made ea polygons III-4 and VII-4. None of 
them has a central symmetry, but III-4 has a full horizontal sym- 
metry, while VII-4 has even no horizontal symmetry. The theoretical 
values of EZ computed for those polygons with C, = 0 are, however, 
much too small, as shown by the crosses on Figure 2. If, however, we 
compute the value of E with C, = 1.9, as if they had central sym- 
metry, the agreement becomes very good. The reason for this is that 
both polygons while not exactly centrally symmetric, are approxi- 
mately so. The deviations from central symmetry fall apparently for 
many individuals within the limits of discriminal error. 

The theoretical values for polygons XII-2, XII-3, and XII-4, as 
shown on Figure 2, have been corrected for end-effects of the rank 
order scale. 

Now consider the polygon of Figure la, and derive from it a set 
of polygons in which the eight sides of the sharpest four angles have 
the same length 1 , while the remaining eight sides have all a different 
length /. All Ohoiee and symmetry relations are the same for all poly- 
gons of each set. By the same argument as before, we would expect 
now the aesthetic value E to be a linear function of > nN, l,, where 


the subscript refers to the L-center. But ~ Ny, ez is ore proportional 
to 81,2 + 81,2. Hence, denoting by A, and N two constants: 
E=A, + Az(l? + L?). (9) 
Figure 3 shows to what extent relation (9) is satisfied. The nega- 
tive slope of the line indicates a large b,, since A, is proportional to 
de (hoe Ope | 
The authors are indebted to Miss Marianne Yampolsky for help 
in recording the experiments and in carrying out the computations. 


This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of the University of Chicago. 
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2. Papers published in the Bulletin cover physicomathematical 
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“1940b”, etc. 

In writing equations, slanted lines should be used wherever pos- 
sible. 
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Drawings should be prepared in a professional manner on white 
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